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Abstract 

We present basic stellar parameters of 99 late-type G giants based on high res- 
olution spectra obtained by the High Dispersion Spectrograph attached to Subaru 
Telescope. These stars are targets of a Doppler survey program searching for extra- 
solar planets among evolved stars, with a metallicity of —0.8 < [Fe/H] < +0.2. We also 
derived their abundances of 15 chemical elements, including four a-elements (Mg, Si, 
Ca, Ti), three odd-Z light elements (Al, K, Sc), four iron peak elements (V, Cr, Fe, 
Ni), and four neutron-capture elements (Y, Ba, La, Eu). Kinematic properties reveal 
that most of the program stars belong to the thin disk. 

Key words: stars: abundances — stars: atmospheres — stars: fundamental 
parameters — stars: late-type 
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1. Introduction 



Since the first extra-solar planet around main sequence star, 51 Peg b was discovered in 
1995 (Mayor & Queloz 1995), more than 500 planets have been revealed in solar neighbourhood 
, providing clues on planet formation process around different types of stars. However, only 
a small fraction of known planet-host stars are intermediate-mass (1.5 ~ 5M ). Properties of 
planets around these stars are of great importance to constrain planet formation theories (e.g. 
Ida & Lin 2004). One way to search for extra-solar planets around intermediate-mass stars is 
monitoring radial velocities of late-type G, K giants. Their spectra have many sharp absorp- 
tion lines, and therefore suitable for getting precise radial velocity data by modern Doppler 
techniques (e.g. Butler et al. 1996). 

The East Asian Planets Search Network (EAPSNet, see Izumiura 2005) which was 
started in 2005 has established an international consortium among Japanese, Korean and 
Chinese researchers. More than 500 late-type G, K giants are being monitored with 2-m class 
telescopes among three countries. Stellar parameters and detailed elemental abundances of 
some of the targets have been analysed by Tekeda et al. (2008) and Liu et al. (2010). As an 
extension of EAPSNet, the Subaru Planets Search Program started monitoring radial velocity 
variations of about 300 late-type G, K giants in 2006. By the end of June 2010, two planets 
(HD 145457b and HD 180314b) have been discovered in this program (Sato et al. 2010). While 
basic parameters and chemical abundance analysis have not been conducted for this sample. 

In this paper we present the first results of stellar parameters and chemical abundances 
for 99 giants based on high resolution spectra obtained with Subaru Telescope. In Section 2, 
we describe the observation and data reduction method. Section 3 gives the derived stellar pa- 
rameters, including the comparison of different methods. Chemical abundances and kinematic 
properties are presented in Section 4 and 5, respectively. In the last section, conclusions are 
presented. 

2. Observation and data reduction 

Targets of Subaru Planet Search Program were selected from Hipparcos Catalogue (ESA 
1997) by the following criteria, 

1. a colour index of 0.6 < B — V < 1.0 and an absolute magnitude of —3 < My ^ 2.5, corre- 
sponding to the range of evolved, late-type G giants. 

2. a magnitude of 6.5 < V < 7.0, which enables follow-up observations using 2m-class tele- 
scopes. 

3. a declination of 5 > —25° so as to be observable in Japan, Korea and China. 

First, each target was observed with the High Dispersion Spectrograph (HDS; Noguchi et al. 
1 see http://exoplanets.org/ 
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2002) equipped to Subaru Telescope three times with an interval of about 1.5 months, for the 
purpose of quickly identifying stars with large radial velocity variations. Then the most likely 
candidates with sub-stellar companions were repeatedly observed with the 1.88m telescope 
at Okayama Astrophysical Observatory (OAO, Japan), 1.8m telescope at Bohyunsan Optical 
Astronomy Observatory (BOAO, Korea), and 2.16m telescope at Xinglong Station (National 
Astronomical Observatories, China). 

HDS has several setups of gratings, and each has different spectrum format and wave- 
length coverage. In the first three runs of this program in 2006, we used StdI2b twice and 
StdI2a once. The spectra taken with StdI2a setup with a wavelength coverage of 4900~7600 A 
are suitable for chemical abundance analysis. Here we present the first results of these 99 stars 
whose spectra were obtained in July, 2006, with a resolving power of R ~ 60,000 and a typi- 
cal signal-to-noise ratio (SNR) of 150-230/pixel. The bias is determined and subtracted using 
the over-scan region on CCD, with an IRAF script provided by HDS web page 2 . Then, the 
reduction of data follows the standard routines. The MIDAS/ECHELLE package is used for 
order definition, flat-fielding, background subtraction, 1-D spectra extraction, and wavelength 
calibration. The radial velocity shift is corrected by fitting profiles of about 80 pre-selected 
lines with intermediate strength. Although the Doppler shift correction is made on the un- 
normalized spectra, the spectra lines are narrow (typically ~ 0.3 A) so that the continumm 
around such a line could be seen as a constant, and nearly do not affect the measured central 
wavelength. The uncertainties of radial velocity values are about 0.3 ~ 0.4 kms -1 using this 
method (see table 8). After that, the one dimension spectra are normalized by a continuum 
fit by a cubic spline function with a smooth parameter. The spline interpolation is determined 
by a set of continuum windows selected from a typical giant spectrum, usually has a density of 
10-15 points per order. 

For the purpose of precise radial velocity measurement, the spectra were taken with an 
iodine vapor cell inserted into the light path of the spectrograph, hence the stellar spectral lines 
located in the wavelength of 5000~6300 A region were heavily mixed with absorption features 
in the I2 spectrum. We only use the red part of the spectra with wavelength longer than 6350 A 
to analyse the chemical abundances. Equivalent widths are measured by fitting the line profiles 
with a Gaussian function for the weak lines. While in strong unblended lines, broad damping 
wings contribute significantly to the equivalent width. The direct integration is used for such 
kind of lines alternatively. 

2 see http://www.naoj.org/Observing/Instruments/HDS/ 
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Fig. 1. Comparison of effective temperatures derived from different colour indices, (a): T c s derived from 
B — V versus that from b — y. (b): T c g derived from V — K versus that from b — y. 

3. Stellar parameters 

3.1. Effective temperature 

Based on the empirical calibration relations given by Alonso et al. (1999, 2001), the 
effective temperatures (T e s) of our program stars are determined using the B — V, V — K, 
and b — y photometric data. The B — V colour indices in Johnson system are obtained by 
B — V = 0.850(i?T — Vt), as described by Perryman et al. (1997), where By and Vr are Tycho 
magnitudes in Hipparcos Catalogue (ESA 1997). The K magnitudes in V — K are obtained by 
converting the K s magnitude in Two Micron All Sky Survey (2MASS) to K magnitude using 
the calibration relation given by Ramirez & Melendez (2004), and V magnitudes are taken from 
Hipparcos Catalogue. There are 32 stars of our sample also have b — y data in catalogue given 
by Hauck & Mermilliod (1998). As a result, we derive the T e g from b — y besides V — K and 
B — V, and compare the results with those from B — V and V — K in figure 1. 

The color excess E(B — V)a is calculated according to Schlegel et al. (1998), with a 
slightly revision described by Arce & Goodman (1999) for those values larger than 0.15 mag. 
Then, the E(B — V) value for each star is calculated as 

E(B-V) = [l-exp(-\Dsmb\/125)]E(B-V) A 

where D is the distance of the star and b is the galactic latitude. Finally, we adopt E(V — K) = 
2.9A8E(B-V) as colour excess for V-K (Schlegel et al. 1998), and E(b-y) = 0.741^(5 -V) 
for b — y (Crawford et al. 1976). 

The mean difference (T c g(B — V) — T c g(b — y)) is 16±56K, while the mean difference 
(T eS (V - K) - T eS (b - y)} is 24±100K. According to Alonso et al. (1999), the uncertainty 
of the effective temperature of giants derived from V — K is estimated to be ±40 K, which is 
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more accurate than that is derived from either B — V (±96 K) or b — y (±70 K). However, most 
stars in our sample are too bright in K s band (< 4 m .5) to get accurate photometry, therefore 
the effective temperatures derived from V — K may not be reliable. Table 1 list the effective 
temperatures derived from three different colour indices for the 32 stars. 



Table 1. Comparison of effective temperatures derived from three different colour indices. 
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Another method is the so-called "excitation equilibrium method" , which forces the de- 
rived iron abundances given by different Fel lines to be independent from their excitation po- 
tentials of lower states (xiow)- After adjusting effective temperatures, the slopes on logA— xi ow 
diagrams for all stars in our sample are smaller than 0.002 dex/eV, and the mean difference 
(T c ff (eq) — T c fj (B — V)) is 44 ± 117K. The effective temperatures derived from this method and 
those from B — V are compared in figure 2. 

3.2. Surface gravity 

Surface gravity (logg) is determined by 

log£ = log£ ±log(^n + 41og(^-J ±O.4(M bol -M boL0 ) (1) 

where M is the stellar mass, and M bo i is the bolometric magnitude defined as 

M bol = V mag + BC - 5 \ogd + 5 - A Y (2) 

where V m3ug , BC, d, and Ay represent the apparent magnitude, bolometric correction, distance, 
and interstellar extinction, respectively. 

The bolometric corrections (BC) are calculated based on estimated effective tempera- 
tures and metallicties, as given in Alonso et al. (1999). The stellar mass is estimated by an 
interpolation of Yale- Yonsei stellar evolution tracks (Demarque et al. 2004, Yi et al. 2003) which 
is based on a new convective core overshoot scheme. Yale- Yonsei evolution tracks can be used 

5 




4400 4800 5200 5600 6000 



T eS (B-V) (K) 

Fig. 2. Comparison of effective temperatures derived from B — V and excitation equilibrium methods. 

for stars from the stage of the pre-main-sequence to the helium-core flash. For comparison, 
we also estimated the masses using the evolution tracks of Girardi et al. (2000), in which the 
helium-core ingnition phase is included. By assuming all the stars are corresponding to the 
post-RGB phase, we find the difference between the two sets of stellar mass is less than O.3M . 
While about half of our program stars present an uncertainty of parallax up to 15%. According 
to equation (1) and (2), this will cause the error of ~ 0.13 dex in logg, and about a factor of 1.3 
in stellar mass. The interstellar extinctions are adopted by A v = 3.1E(B — V), and Hipparcos 
parallaxes are used to determine the absolute magnitudes (My). 

We also adopt another method which is independent from parallax data to determine 
the surface gravities. As the Fell lines are more sensitive to surface gravities in cool stars than 
Fel lines, log g can be determined by forcing the abundances of Fell lines to be the same value 
as those given by Fel lines. Despite several uncertainties such as non-local thermodynamic 
equilibrium (NLTE) effects and very few Fe II lines, it is still meaningful to compare different 
sets of logg values. As shown in figure 3, the scatter tends to increase with decreasing surface 
gravities. 

3.3. Atomic data 

The log gf values for selected spectral lines are taken from various literatures, and listed 
in table 2. For Fe I lines, we use Oxford oscillator strength data (Blackwell 1982a,b, 1986). For 
the four heavy elements (Y, Ba, La, Eu), the logg/ values are taken from Hannaford (1982), 
Weise & Martin (1980), Luck & Bond (1992), and Biemont et al. (1982), respectively. For 
other elements, we use the logg/ values given by Chen et al. (2000). Although the atomic data 
are compiled from different references, it does not affect the results systematically. Moreover, 
the same program and values of oscillator strength are also adopted by measuring the EWs in 
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Fig. 3. Comparison of surface gravities obtained by Hipparcos parallaxes (Y axis) and ionization balance 
of Fel and Fell lines (X axis). 

Solar Atlas (Kurucz et al. 1984) to determine the solar abundances, so the final results are the 
differential values relative to the Sun. 

Table 2. Atomic line data. 



Aair (A) 


Element 


Xlow 


log.9/ 


E.~, 


6327.604 


Ni I 


1.68 


-3.150 


2.5 


6335.337 


Fe I 


2.20 


-2.177 


1.2 


6336.830 


Fe I 


3.69 


-0.856 


1.4 


6344.155 


Fe I 


2.43 


-2.923 


1.3 



3.4- Metallicity and microturbulence velocity 

For most of our program stars, the initial metallicity values are set to [Fe/H]=0.0. For 
those stars that have been studied, the initial values are taken from previous literatures (da 
Silva et al. 2006, Mishenina et al. 2006, Luck & Heiter 2007, and Tekeda et al. 2008). The 
abundances of chemical elements are determined based on the model atmospheres interpolated 
by a plane-parallel, homogeneous and local thermodynamic equilibrium (LTE) model grid by 
Kurucz (1993). Chemical abundances are calculated with ABONTEST8 program supplied 
by Dr. P. Magain (Liege, Belgium). It calculates the theoretical equivalent widths from the 
atmospheric model, and matches them with observed values. Several broadening mechanism 
have been taken into account, including natural broadening, thermal broadening, van der Waals 
damping broadening, and microturbulent broadening. The microturbulence £ t is determined 
by forcing the iron abundance values given by different Fe I lines to be independent from their 
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EWs and hence with zero slope on log A-EW diagram. Only those Fe I lines with 10 mA 
< EW < 110 mA are used. 

3. 5. Summary 

The derived physical parameters of our program stars are summarized in table 3 and 4. 
In table 3, effective temperatures are derived from B — V and and surface gravities are derived 
from Hipparcos parallaxes. While in table 4, effective temperatures and surface gravities are 
derived from excitation equilibrium and and ionization balance, respectively. Although two 
different sets of stellar parameters are given, we finally adopt those listed in table 3 in the 
following analysis. 



Table 3. Physical parameters of program stars. Effective temperatures are derived from B — V and calibration relation of 
Alonso et al. (1999). Surface gravities are derived from Hipparcos parallaxes. 
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-0.351 
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Table 4. Atmospheric parameters of program stars. Effective temperatures are 


derived from excitation equilibrium method. 




Surface 


gravities are derived from ionization balance method. 
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In figure 4 log(L/L Q ) versus T e g is plotted for our program stars. It is shown that most of 
our samples have luminosities 1.5 < logL/L© < 2.0, and effective temperatures 3.67 < logT/T c g < 
3.70. In figure 5, the derived atmospheric parameters of our stars (filled circles) are compared 
with those derived by Tekeda et al. (2008) (open circles) and Liu et al. (2010) (crosses). The 
effective temperature {T e g) tends to increase with increasing surface gravity (log g), and the 
microturbulent velocity (£ t ) tends to decrease with higher surface gravity (logg). 

3.6. Consistency check with literatures 

In order to check the consistency of our atmospheric parameters with previous studies, 
the results are compared with those given by da Silva et al. (2006), Mishenina et al. (2006), 
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Fig. 4. Luminosity versus effective temperature of our program stars, together with the Yale-Yonsci 
evolution tracks (Dcmarque et al. 2004) of stars with different masses (1M Q , 2Af Q , 3M Q , and 5M 0) ) and 
metallicity (z=0.04, 0.02, and 0.01, corresponding to [Fe/H] =+0.3, 0.0, and -0.3, respectively). 
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Fig. 5. Relations of atmospheric parameters (log<? vs. T e g, [Fe/H] vs. T e g, £t vs. logg) derived from 
B — V. Filled circles: this study, open circles: Tekeda et al. (2008), and crosses: Liu et al. (2010). 
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Hekker & Melndez (2007), Luck & Heiter (2007), Tekeda et al. (2008), and Sato et al. (2010). All 
literatures above have determined the parameters of giants based on high resolution spectra. 
There are totally 7 common stars in our sample with those literatures, as listed in table 5. 
Figure 6 compares the effective temperatures, surface gravities, and metallicity of these stars. 
Our derived temperatures are well consistent with others, with AT er r ~4±65K lower than their 
values. While the surface gravity difference Alogg is ~ 0.14 ±0.33 dex lower compared with 
others' results, and the metallicity difference A[Fe/H] is ~ 0.09 ±0.07 dex lower than those 
from literatures. 

Table 5. Common stars in this study and other literatures 



HD 


Literature 


145457 


Sato et al. 2010 


161502 


Luck & Heiter 2007 


179799 


da Silva et al. 2006 


180314 


Sato et al. 2010 


185351 


Tekeda et al. 2008 




Hekker & Melndez 2007 


188993 


Luck & Heiter 2007 


196134 


Mishenina et al. 2006 



4. Chemical abundances 

We use the atmospheric parameters listed in table 3 to derive the [X/Fe] ratios of 14 
elements (Al, Ba, Ca, Cr, Eu, K, La, Mg, Ni, Sc, Si, Ti, V and Y). Those trends are plotted 
in figure 9 to figure 13 (filled circles), together with the trends from Takeda et al. (2008, open 
circles) as a comparison. The measured equivalent widths of the selected lines of our program 
stars are listed in table 9, which is only available in the electronic form. 

4-1. Error analysis 

The uncertainties of chemical abundances are estimated by changing the atmospheric 
parameters in a reasonable range. Table 6 shows the abundances differences due to deviations 
of the effective temperature of 100 K, the surface gravity of 0.15 dex, the iron abundance of 
0.1 dex, and the microturbulent velocity of 0.1 kms -1 . For most of the chemical elements, the 
uncertainties are less than 0.1 dex, except of Ti and V. The uncertainties of these two elements 
can be as high as 0.15 dex. 

Another source of uncertainties is the error of equivalent widths caused by mixture of the 
intrinsic stellar lines and the I2 lines, since the spectra of our program stars are 12-superposed 
for precise radial velocity measurements. To check the uncertainties caused by iodine lines, we 
use the pure stellar spectra of HD 145457 taken in OAO to measure the equivalent widths, and 
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logg (this work) 

Fig. 6. Comparison of atmospheric parameters of common stars in this study and other literatures. 
Upper panel: effective temperature, middle panel: metallicity, bottom panel: surface gravity. 

determine the chemical abundances. 

In Doppler surveys, spectra of fast-rotation B type stars are often obtained with I2 cell 
to deconvolute instrumental profiles (IPs, e.g. Butler et al. 1996). The spectra of such stars are 
featureless, and thus provide a simple way to check how the stellar spectra are contaminated 
by I2 absorption lines at different wavelengths in this study. The ^-superposed spetra of a B 
star (HR 5685, B8V), I 2 -superposed spetra of HD 145457, and the pure spectra of HD 145457 
around 6400 A and 6700 A are plotted in figure 7. Around 6400 A, the integrated equivalent 
width of I2 lines is about 6mA/A, which causes the uncertainty of equivalent width of about 
3 ~ 4 mA for a typical spectral line in our ^-superposed spetrum in this region. While around 
6700 A, the integrated equivalent width of I2 lines is about 4mA/A, 30% smaller than that 
around 6400 A. The uncertainty of equivalent width caused by I2 lines is estimated to be 2mA. 
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Table 6. Estimated errors for abundances of a typical late-type giant HD109305, with T e ff = 4663 K, logg = 2.44, 
[Fe/H] = -0.30, and ft = 1.77kms _1 , taken from table 3. 
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Fig. 7. Comparison of the 12-superposed spectrum (taken in Subaru) and the pure spectrum (taken in 
OAO) for HD 145457 around 6400A and 6700 A. The calculated RMS residuals are 1.2% around 6400A 
and 0.7% around 6700 A. Two portions of 12-superposed spectrum for a fast rotation, B-type star (HR 
5685, taken in Subaru) are also plotted in the upper panel. 
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Fig. 8. Comparison of equivalent widths (EWs) measured using ^-superposed spectra taken in Subaru (X 
axis) and pure spectra taken in OAO (Y axis) . The spectral lines are separated into three groups according 
to their wavelengths. The filled circles are 12 lines with wavelengths A < 6500 A. The open circles are 44 
lines with wavelengths 6500 A < A < 7000 A. The crosses are 31 lines with wavelengths A > 7000 A. The mean 
differences (EW U - EW puIe ) are -0.05 ± 5.22mA, -0.33 ± 4.33mA, and +1.48 ± 6.38mA, respectively. 

Furthermore, we measure the equivalent widths using the ^-superposed spctrum and 
pure spectrum of HD 145457. The spectral lines are separated into three groups, with wave- 
lengths A < 6500 A, 6500 A < A < 7000 A, and A > 7000 A, respectively. The results are com- 
pared in figure 8. We find that the mean differences (EWu — EW pure ) of three groups are 
-0.05 ± 5.22 mA, -0.33 ± 4.33 mA, and +1.48 ± 6.38 mA, which infers our EW values are not 
significantly influenced by I 2 absorption lines. The chemical abundaces determined by pure 
stellar spectra are compared with our results using Subaru spectra in table 7. It is shown 
that the uncertainties of chemical abundances caused by mixture of I 2 lines are range from 
0.02 ~ 0.09 dex, which are smaller than the uncertainties caused by errors of atmospheric pa- 
rameters. Therefore, we concluded that the weak I 2 absorption lines in the red part of our 
spectra nearly do not affect the chemical abundances measurements. 

4-2. a-elements (Mg, Si, Ca, Ti) 

The a elements are mainly produced in SNe II nucleosynthesis (Woosley & Weaver 1995) 
and show enhancement in metal-poor stars. In figure 9, we plot [X/Fe] ratios of Mg, Si, Ca, and 
Ti I versus [Fe/H], together with the trends take from Tekeda et al. (2008). Except for the trend 
of [Til/Fe] versus [Fe/H] with a large scatter, three of the other a elements exhibit enhancement 
towards lower metallicity. However, their decreasing trends with increasing metallicity are 
slightly different. Both the [Mg/Fe] and [Ca/Fe] show a turn off trend towards a flat pattern 
at [Fe/H] ~ —0.2, which is in good agreement with previous studies on both giants (Tekeda et 
al. 2008, Liu et al. 2007) and dwarfs (Chen et al. 2000), while the [Mg/Fe] turns to decrease at 
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Table 7. Comparison of chemical abundances of HD 145457 determined by pure stellar specta taken in OAO, and I2 
superposed spectra taken in Subaru. 



AIX/H1 


[X/Hli 2 - fX/Hl Duro 


A[FeI/H] 


0.02 


A [Fell /Hi 


0.04 




0.04 


A[A1/H] 


0.03 


AfSi/Hl 

L / J 


0.04 


A[Ca/H] 


0.09 


A[Sc/H] 


0.06 


AfTi/Hl 

L / J 


0.04 


AfV/Hl 


0.02 


AfCr/Hl 


0.07 


A[Ni/H] 


0.07 


A[Y/H] 


0.02 


A[Ba/H] 


0.02 


A[La/H] 


0.03 


A[Eu/H] 


0.02 


A[K/H] 





[Fe/H] > 0. Despite an increasing trend with decreasing metallicity for [Fe/H] < —0.4, [Si/Fe] 
shows a larger scatter for higher metallicity, leaving the flat pattern shown in Tekeda et al. 
(2008) (open circles in figure 9) unclear. Oxygen is another important a element for giants. 
But in our spectra, the [O I] forbidden lines at 6300 A are not covered by the CCD, and the 
lines at 6363 A are very weak and closed to the edge. To make sure our results are reliable, the 
oxygen abundance is not included in this study. 

4-3. Odd-Z light elements (Al, K, Sc) 

As shown in figure 10, [Al/Fe] shows a decreasing trend with increasing metallicity for 
[Fe/H] < —0.2, and a solar pattern towards higher metallicity. Our results are consistent with 
Liu et al. (2007). Chen et al. (2000) found a rather steep upturn of [Al/Fe] at [Fe/H] = -0.2 
in a sample of dwarfs. However, in our results the upturn is rather smooth. [K/Fe] exhibits a 
larger dispersion, because only very few number of lines are available. The trend of [Sc/Fe] is 
similar with that of [Ti/Fe], with both of them showing a flat pattern when data with larger 
dispersion are excluded. This result is also consistent with Tekeda et al. (2008). 

4-4- Iron-peak elements (V, Cr, Ni) 

The iron-peak elements are believed to be mainly produced in SNe la, and should have 
the same abundance patterns with iron. In figure 11, we plot [X/Fe] versus [Fe/H] for three iron- 
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Fig. 9. Abundance ratio [X/Fe] for four a elements (Mg, Si, Ca, and Ti I) against [Fe/H] (filled circles), 
together with the results from Takeda et al. (2008, open circles). 

peak elements (V, Cr, and Ni). [V/Fe] shows an increasing trend with increasing metallicity, 
which is contrast to Takeda's results (plotted as open circles in the same figure). The hyper fine 
structure (HFS) effect has been considered in our work. The adopted HFS data is taken from 
Kurucz 3 . Liu et al. (2007) suggests the HFS effect can lead to a correction as large as 0.5 dex 
for vanadium abundance in giants. However, we find that the corrections due to introduction 
of HFS effect are between -0.002 dex and 0.012 dex, and can be neglected (see figure 12). The 
discrepancy is due to the difference of HFS effect from line to line, and the vandium lines we 
used in this work are different from those in Liu et al. (2007). 

The [Cr/Fe] trend in our sample is consistent with Takeda's, except for a few stars with 
higher [Cr/Fe] values. For these stars, the relation of [Ni/Fe] vs. [Fe/H] keeps the solar pattern 
in the range of —0.8 < [Fe/H] < —0.1, while shows overabundance for higher metallicity and 
hence systematically higher than Takeda's result. The possible upturn at [Fe/H] ~ 0.0 has also 
been found in giants (e.g. Liu et al. (2007)) and dwarfs (e.g. Edvardsson et al. 1993). 

3 http:/ /kurucz. harvard.edu/linelists. html 
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Fig. 10. Abundance trends of [X/Fe] against [Fe/H] for three odd-z elements (Al, K, and Sc, filled circles). 
Results of [Sc/Fe] from Tekeda et al. (2008) are plotted as open circles. 




Fig. 11. Abundance trends of [X/Fe] against [Fc/H] for three iron-peak elements (V, Cr, and Ni, filled 
circles). Results from Tekeda et al. (2008) are plotted as open circles. 
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Fig. 12. The difference in vandium abundances obtained by neglecting HFS minus those by considering 
HFS versus [Fe/H]. A[V/Fe] on Y axis represents logV non -HFS — log'V HFS . 
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Fig. 13. Abundance trends of [X/Fe] against [Fe/H] for four heavy elements (Y, Ba, La, and Eu). (filled 
circles: this work; open circles: Tekeda et al. 2008's work) 

4-5. Heavy elements (Y, Ba, La, Eu) 

Those elements heavier than iron are mainly produced via the neutron capture process. 
Two major mechanisms are generally considered: the r-process and s-process, corresponding to 
neutron rich and neutron poor environment, respectively. In our analysis, one r-process element 
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Fig. 14. [La/Eu], [La/Y], [La/Ba] as a function of [Fe/H] (from top to bottom). 

(Eu), one light s-process element (Y) and two heavy s- process elements (Ba, La) are included. 
Their abundances trend against [Fe/H] are plotted in figure 13. 

Despite a large scatter as often observed in field stars, the ratio of [La/Eu], [La/Y], and 
[La/Ba] may provide important information of the Galactic chemical evolution history. As a 
typical r-process element, europium is mainly produced in core-collapse supernovae (SNe) with 
masses 8M Q < M < 10 M and the pattern still remains uncertain. Previous work (Simmerer et 
al. 2004) of La and Eu revealed the s-process may be active as early as [Fe/H]=-2.6. In figure 
14, [La/Eu], [La/Y], and [La/Ba] are plotted as functions of [Fe/H] (from top to bottom). From 
our plot, a decline of [La/Eu], although not very clearly, can be seen in lower metallicity down 
to -0.7. In the convection zone of low mass AGB stars where the slow neutron capture reactions 
take place, the neutron flux per seed nuclei is inversely proportional to the density of iron-peak 
nucleus. The metal-poor stars produce light s-process elements (such as Y) more efficiently 
than the heavy s-process elements (such as Ba, La). The [La/Y] versus [Fe/H] in the middle 
panel of figure 14 shows an marginal increasing trend with decreasing metallicity, consistent 
with the theoretical prediction (e.g. Busso et al. 1999). [La/Ba] shows a flat trend with [Fe/H] 
in the bottom panel of figure 14. This is consistent with the theoretical prediction that they 
are produced through the same producing mechanism. 

5. Kinetics Parameters 

We derived the kinetics parameters (J/lsr, ^lsr; Wlsr) based on the method given 
by Johnson & Soderblom (1987). The radial velocities (RVh e u ) are measured based on our 
spectra, and corrected to heliocentric velocities. Positions, parallaxes, and proper motions are 
taken from Hipparcos data. By adopting a solar motion of (U,V,W) Q = (— 10. 00 ±0.36, +5.25 ± 
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Fig. 15. Toomre diagram of program stars. Thin disk stars (TD/D < 1.0) and thick disk stars 
(TD/D > 1.0) arc plotted as filled and open circles, respectively (See discussion in Sect. 5). Dotted 
lines indicate constant total velocities |«|lsr = v^lsr + ^lsr + ^lsr m s ^ e P s °f 50 kms -1 . 

0.62, +7.17 ± 0.38) 4 in kms -1 , given by Dehnen & Binney (1998), the Galactic velocity com- 
ponents U, V, and W of our program stars are corrected to local standard of rest (LSR). It 
is hard to clarify whether a star in solar neighborhood belongs to the thin disk or the thick 
disk. We apply the kinematic method proposed by Bensby et al. (2003) to calculate the rela- 
tive probability for the thick-disk-to-thin-disk (TD/D) membership. In Bensby et al. (2003), 
thin disk and thick disk stars are thought to follow Gaussian distributions in Galactic velocity 
space (U, V, W). By assuming different asymmetric velocity drifts and dispersions, the relative 
probability of belonging to different populations can be calculated. For instance, a star with 
TD/D = 10 means it is 10 times more likely to belong to the thick disk than to the thin disk. 
In our sample, one star (HD 115903) with TD/D > 10.0 and two stars (HD 138425 and HD 
176973) with 1.0 < TD/D < 10.0 are more likely to be a thick disk star rather than a thin disk 
star. Those stars are plotted as open circles on Toomre diagram in figure 15, and listed with a 
mark 'TD' in the last column of table 8, as well as their kinetics parameters mentioned above. 

6. Conclusion 

In this work we determine the atmospheric parameters, kinematic properties and chem- 
ical abundances of 99 late-type giants based on high resolution spectra obtained in Subaru 
Planet Search Program, covering the metallicity range —0.8 < [Fe/H] < +0.2. We get two sets 
of stellar parameters derived from different methods. The photospheric chemical abundances 
of 15 elements, including four a-elements (Mg, Si, Ca, Ti), three odd-Z light elements (Al, K, 

4 In this work the Galactic velocity component U is defined to be positive towards Galactic anticenter. 



19 



Table 8. Kinetics parameters of program stars 



HD 


-KV'helio 




T T 

VXSR 


WlSR 


LD/D 




(kms -1 ) 


(kms -1 ) 


(kms- 1 ) 


(kms" 1 ) 




lUUUOO 


ft fin j- n Q7 
0.09 ± U.o 1 


z.y ± i.o 


-10.4 dz Z.D 


lo.O ± U. ( 


U.U1 


101853 


2.81 ± 0.36 


-30.9 ± 3.4 


6.7 ± 0.7 


17.3 ± 1.3 


0.01 


103690 


-19.28 ± 0.35 


-58.6 ± 3.8 


5.0 ± 1.1 


8.9 ± 1.7 


0.02 


105475 


0.35 ± 0.35 


-80.9 ± 11.8 


-0.4 ± 1.3 


16.9 ± 1.7 


0.06 



Sc), four iron peak elements (V, Cr, Fe, Ni), and four neutron-capture elements (Y, Ba, La, 
Eu) are determined. The kinematic parameters are calculated, and most of our sample are thin 
disk stars. From the results of abundances, we conclud that the abundance ratios [Mg/Fe] and 
[Ca/Fe] show an increasing trend with decreasing metallicty for —0.8 < [Fe/H] < —0.2, while 
flatten out towards higher metallicity. A decreasing trend of [Mg/Fe] with increasing metallicity 
is detected for giants at [Fe/H] > 0. The upturn of [Al/Fe] versus [Fe/H] at [Fe/H] = -0.2 in 
giants is not as steep as that in dwarfs (e.g. Chen et al. 2000). The [Ni/Fe] versus [Fe/H] shows 
an upturn at [Fe/H]=0.0 and an overabundance towards higher metallicity. The [La/Y] ratio 
shows an increasing trend with decreasing metallicity, which is consistent with the prediction 
of AGB nucleosynthesis. 
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